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Abstract. We study the evolution of two mixed Dirac neutrinos in an external axial- vector 
field. The dynamics of this system is described in the framework of relativistic wave equations 
approach on the basis of the known solutions of the Dirac equation in an external field with a 
given initial condition. The general solution of the initial condition problem, exactly accounting 
for the external field, is obtained. Then we consider special form of the external field which 
corresponds to the standard model neutrino interactions with the background matter. In this 
limit we derive the transition probability and compare it with the previously known results. 



The problem of neutrino oscillations in various external fields attracts considerable attention 
mainly because the amplification of neutrino oscillations in background matter (MSW effect, see 
Refs. [HE]) is the most plausible explanation of the solar neutrino deficit. We discussed neutrino 
flavor and spin-flavor oscillations in our recent works [U EJ [6] in frames of the relativistic wave 
equations approach (classical field theory) which exactly accounts for external fields. The present 
paper continues our previous efforts. 

Let us discuss the evolution of two mixed flavor neutrinos v\, A = a, j3, interacting with the 
external field axial- vector field fy\i- We assume that the matrix (/vy) is generally not diagonal 
in the indexes A, A'. Note that the similar external fields were recently discussed in Ref. [7] (see 
also references therein). The Lagrangian describing the dynamics of the considered system is 

X=a,f3 \,\'=a,P 



where (m\\>) is the mass matrix of flavor neutrinos and 7^ = 7^(1 + 7 5 )/2. Note that the 
matrices (m\\>) and (fxy) are generally independent. We supply the Lagrangian (TjQ) with the 
initial conditions, v a (r,0) = and ^(r, 0) = £(r), where £(r) is a given function. This situation 
is implemented, e.g., if v a = v^, up = v e and we study oscillations of solar neutrinos, when only 
electron neutrinos are presented initially. 

To examine the evolution of the system (TjQ) we introduce the mass eigenstates ip a (r, t), a = 1,2, 
by means of the matrix transformation, v\ = ^ U\ a ip a . In case of two flavor neutrinos the matrix 
(U\ a ) is the rotation matrix in the flavor basis, U = exp(— io"2#), where 9 is the vacuum mixing 
angle. The effective Lagrangian for the mass eigenstates is obtained from Eq. (pQ), 

CfaM = J2 Mh^-maWa- 9abi>a^b, (2) 

o=l,2 a,fe=l,2 

where m a the mass is fermion tp a and (g^ b ) = U~ 1 (f^y)U is the external axial- vector field 
expressed in the mass eigenstates basis. 



The Dirac equations which result from Eq. ([2]) have the form, 



iip a = H a tp a + Vtpb, a, 6 = 1,2, a / 6, 



(3) 



where 7i a = — ictV + (3m a + p-y^gii, V = /?7^g M , ga = giia and g^ = g± 2 - The general solution to 
Eq. ([3]) can be presented in the following way (see Refs. [UEJE]): 
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where ai and 6a are the undetermined non-operator coefficients, which are generally time 
depedent. The energy spectrum Ea hi case of non-moving and unpolarized mater, which 
corresponds to g a = g 



0, was found in Refs. [HI [§] 
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In Eq. © we introduce the quantity g a =g®. The basis spinors u a ^ and vffi in Eq. Q are the 
eigenvectors of the helicity operator (Sp)/|p|, with the eigenvalues £ = ±1, and can be also 
found in the explicit form in Refs. [HJ[9]. 

Let us assume the initial wave function in the form, £(r) = e lkr £o> where k = (0, 0, k). If we 
study relativistic neutrinos with k 3> mip, the normalized basis spinors in Eq. are 
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One can take that £o = u~ . It is easy to check that (1/2)(1 — £3)^0 = Co- Therefore £(r) 
describes a neutrino propagating along the z-axis, with the spin directed opposite to the particle 
momentum. Note that the subscript a is omitted in Eq. since we neglect small terms 
~ m a /k <C 1. 

With help of the obvious identities, {u~\V\u~) = (v + \V\v + ) = g° (all other matrix elements 
of the potential V vanish), which result from Eq. ([6]), as well as using the technique developed 
in Refs. 0E] we get the ordinary differential equations for the functions a~(t), 



ia„ 



a b # ex P K E a ~ E b )*]> a,6 = l,2, a^b, (7) 

where g = g°. Note that the equations for the functions b£(t) are obtained analogously. The 
solution to Eq. © can be expressed in the form (see, e.g., Refs. E]) 



a -(t) =Fa~[(Q) + Ga,2 (0), 
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exp (iut/2), 



a~{t) =F*a 2 (0) - G*a7/(0), G = - sin exp (iart/2), 



(8) 



where = ^yg^^(uj/2)^ and a; = -E] - — -ET . 

Using the identity (v + (g) Co = [see Eq. as well as Eqs. and (|HJ) we arrive to the 
wave function of the neutrino u a , 

/ x ■ , • i \ • *-\ \a cos 26* + (uj/2) sin 201 . _ fm a \ ,„ s 
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where E = (E^ +£ , 2 ")/2. Note that Eq. ([9]) is the most general one which is valid for the external 
axial- vector fields f£ x , of arbitrary strength. Let us however discuss one of the applications of the 
obtained result. We consider diagonal matrix f^ x , = f^5\y which corresponds to the standard 
model neutrino interactions with the non-moving and unpolarized background matter, i.e. we 
take that f*A = 0. Note that the values of /jj for various channels of neutrino oscillations can 
be found in Ref. [TDj In this case we have g = — sin2#A/°, where A/ = (/° — /«)/2. In 
the low density matter limit, g a <C k, one reads lu/2 ~ <l?(fc) + cos2#A/° [see Eq. (j5J)], where 
= 5m 2 /(4k) and 5m 2 = m\ — to 2 . 
The transition probability can be calculated on the basis of Eq. ([9]) as 



where 



P,.. -,..(/) - |/^C./)|- ~ ■ (.10) 



A = — — _ -, - = a/ ifjfe) cos 20 + A/0] 2 + $2(fe) sin 2( 2 0). (11) 

[$(fc)cos26> + A/°] 2 + $ 2 (£;)sm 2 (260 L V L v ' v ; v ; v ' 

We can observe that Eqs. (|10|) and reproduce the famous formula for the neutrino oscillations 
probability in the background matter (see Refs. [H12])- 

In conclusion we mention that neutrino flavor oscillations in background matter with non- 
standard interactions have been studied in frames of the relativistic wave equations approach 
(classical field theory). Neutrino interactions with matter are known to be equivalent to the 
presence of an external axial-vector field. In the limit of ultrarelativistic neutrinos interacting 
with non-moving and unpolarized matter we have received the most general expression for 
the neutrino wave function [Eq. ([9"])], which exactly takes into account the external axial- 
vector field. Using this result and considering standard model neutrino matter interactions 
we could reproduce the transition probability [Eqs. (|10p and (jlip ] accounting for the resonance 
enhancement of neutrino oscillations - the MSW effect. Note that with help of the method 
used in the present work we have improved the results of our recent paper [5]. In contrast to 
that work, where the case of only low density matter was discussed, now we could obtain the 
transition probability valid for arbitrary matter density. 
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